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Abstract 

I consider the problem of self-oscillatory systems undergoing a homogeneous Hopf bifurcation 
when they are submitted to an external forcing that is periodic in time, at a frequency close 
to the system's natural frequency (1:1 resonance), and whose amplitude is slowly modulated in 
space. Starting from a general, unspecified model and making use of standard multiple scales 
analysis, I show that the close-to-threshold dynamics of such systems is universally governed by a 
generalized, complex Ginzburg-Landau (CGL) equation. The nature of the generalization depends 
on the strength and of other features of forcing: (i) For generic, sufficiently weak forcings the CGL 
equation contains an extra, inhomogeneous term proportional to the complex amplitude of forcing, 
as in the usual 1:1 resonance with spatially uniform forcing; (ii) For stronger perturbations, whose 
amplitude sign alternates across the system, the CGL equation contains a term proportional to 
the complex conjugate of the oscillations envelope, like in the classical 2:1 resonance, responsible 
for the emergence of phase bistability and of phase bistable patterns in the system. Finally I show 
that case (ii) is retrieved from case (i) in the appropriate limit so that the latter can be regarded 
as the generic model for the close-to-threshold dynamics of the type of systems considered here. 
The kind of forcing studied in this work thus represents an alternative to the classical parametric 
forcing at twice the natural frequency of oscillations and opens the way to new forms of pattern 
formation control in self-oscillatory systems, what is especially relevant in the case of systems that 
are quite insensitive to parametric forcing, such as lasers and other nonlinear optical cavities. 
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I. INTRODUCTION 



The temporal periodic forcing of spatially extended, self-oscillatory systems is a clas- 
sical method to control and excite the formation of spatial patterns in such systems. 
This kind of forcing admits a universal description when the system is operated near 
the oscillation threshold and forcing acts on an n : m resonance, defined by the relation 
cjf = (n/m) (uq + Sou) between the external forcing frequency Uf and the natural frequency 
of oscillations uq, where n/m is an irreducible integer fraction and Su is a small mistuning. 
In such case oscillations settle down in the system so that any of its variables takes the form 
Re [ku (x, t) e J (^o+' 5w )'J ; to the leading order, where n is a complex constant and the slowly 
varying complex amplitude of the oscillations u verifies the following generalized complex 
Ginzburg-Landau equation (CGLE) jl, 2] 

d t u = a\u + a 2 V 2 w + a 3 |m| 2 u + a^u" 1 * 1 , (1) 

where u stands for the complex conjugate of u, ^1=1,2,3 are complex coefficients and 04 is 
proportional to the m-th power of the forcing amplitude. Equation (JT|) is valid in principle 
for perfectly periodic forcings and has been introduced making use of elegant symmetry 
arguments [l| as well as has been derived in specific contexts, such as in chemistry 3], 
making use of standard multiple scales analysis 4j, . 

As commented the validity of ([1]) requires that the system is close to the Hopf bifurcation 
and, additionally that the amplitude of the external forcing be either constant or slowly 
varying in time and/or space, see {5] for the case n = 1. Here I show that Eq. ([I]) with 
n = 2 applies as well to systems forced at the 1:1 resonance (hence at n = 1) when the 
forcing amplitude varies on space on a long (but not too long) scale, to be defined formally 
below, and the spatial variation involves a sign alternation of the amplitude. This implies 
that the system becomes phase bistable, which is one of the salient features of Eq. (CQ) with 
n = 2 ??, in contrast to the phase locking (to just one phase value) predicted in the case 
n — 1. The intuitive picture is that, as the system is resonantly forced at a 1:1 resonance the 
phase of the oscillations tend to lock to a single value that depends in particular on the phase 
of the forcing. In our case however forcing displays two opposite phases (two signs) that 
are alternating across the system. The system's oscillations have then two reference phases 
and can lock to two dynamically equivalent values, leading to phase bistability. Clearly this 
must require that the typical spatial scale over which forcing varies is short as compared 



with the typical spatial scale of the unforced system as otherwise the system could lock 
locally to the forcing phase at that region. The kind of forcing considered in this work gets 
inspiration from another 1:1 resonant forcing, called rocking [6|, in which the amplitude of 
forcing is uniform in space but its sign alternates along time. Rocking has been considered 
both theoretically and experimentally in several systems [7- loj|. 

The analysis presented here is based on the technique of multiple scales and generalizes 
the concept of resonant forcing as it considers periodic forcings in time that are nonuniform 
across the system on an apparently nonresonant scale. The derivation allows a rigorous 
simplified study of both spatially periodic forcings as well as spatially noisy forcings within 
a 1:1 resonance, whenever they verify the conditions imposed below. Generalizations to 
other resonances (n / 1 or m / 1), although cumbersome, can be made straightforwardly 
following the lines of the present derivation. 

The kind of forcing studied here thus represents an alternative to the classical parametric 
(or 2:1 resonant) forcing, what is especially relevant in systems that are insensitive to the 
latter, like most nonlinear optical systems. In particular this allows the emergence of phase 
bistability in the system (as the term a^u breaks the usual continuous phase symmetry of 
the classical CGLE down to the discrete one u — > —u [l|) and of phase bistable patterns, 
like phase domains and phase domain walls, rolls, hexagons and localized structures. 



II. MODEL AND NOTATION 



We consider a generic one-dimensional system described by N real dynamical variables 
{Ui (x, t)} i=1 whose time evolution is governed by the following set of real equations written 
in vector form, 

dtU (x, t) = f (ji, a;V) + V (ji, a) ■ V 2 U, (2) 

where f is a sufficiently differentiate function of its arguments, V 2 = <9 2 is the one- 
dimensional Laplacian and T> (/1, a) is a diffusion matrix. This is the simplest dependence 
on derivatives in space-translation invariant systems at the time it corresponds to actual 
systems of most relevance like, e.g., reaction-diffusion and nonlinear optical systems. \x is 
the bifurcation parameter and a(x,t) is the forcing parameter, which is allowed to vary 
on time and space. Physically a may represent either an independent parameter, or the 
modulated part of any other parameter. 
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We assume that in the absence of forcing (a = 0) Eq. (J2J) supports a steady, spatially 
homogeneous state U = U s (/x) (d t U s = d x XJ s = 0), which looses stability at /i = /i giving 
rise to a self-oscillatory, spatially homogeneous state. In other words, we assume that the 
reference state U s suffers a homogeneous Hopf bifurcation at /i = /io- We wish to study 
the small amplitude oscillations that form in the system close to the bifurcation when the 
arbitrary parameter a is periodically modulated in time with a frequency close to that of 
the free oscillations and is modulated in space on a long spatial scale to be defined below. 

For the sake of convenience we introduce a new vector 

u(r,t) = U(M)-U s , (3) 

which measures the deviation of the system from the reference state, in terms of which we 
rewrite Eq. (T2J) as a Taylor series, 

<9 t u = F(fi,a) + J(fi,a)-u + V (//, a) ■ V 2 u 

+ K (/j, a; u, u) + L (//, a; u, u, u) + h.o.t., (4) 

where h.o.t. denotes terms of higher order than 3 in u. As usual, and as we show below, these 
h.o.t. have no influence near the bifurcation whenever it is suprecritical, what we assume. 
The different elements of the expansion (T4J) are defined as 

F 0, a) = f s , 
Jij (p, ") = [dfi/dUj^ , 
K(/i,a;a,b) = |E5=i [d^/dUidUj]^^, 
L (/i, a; a, b, c) = ^J2?j,k=i [^ 3f '/dU i dU j dUk] a aib j c k , 

where a, b, c are arbitrary vectors and the subscript "s" denotes U = U s (/i). Vector F (/i, a) 
is subjected to the condition 

F(>,0) = 0, (5) 

since in the absence of forcing (a = 0) the reference state u = is a steady state of Eq. (T4j) 
by hypothesis. J is a matrix and vector K (L) is a symmetric and bilinear (symmetric and 
trilinear) form of its two (three) last arguments. 
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III. THE HOPF BIFURCATION OF THE UNFORCED SYSTEM 



In the absence of forcing (a = 0) the stability of the reference state against small pertur- 
bations 5u is governed by the following equation 

d t 5u = J(ji,0)-6u + V (ji, 0) ■ V 2 5u, (6) 

obtained upon linearizing Eq. (j4]) for a = with respect to 5u. The general solution to Eq. 
([H]) is a superposition of plane waves of the form 

5u (r, t) = J2j w j ex P (A^) exp (ikjx) , (7) 

with 

AjW^Mfakfi-w,, (8) 

M (ji, k 2 ) = J (/i, 0) - k 2 V (fj,, 0) , (9) 

hence eigenvalues and eigenvectors of matrix Ai depend on k only through k 2 as Ai does 
(this is a consequence of the assumed spatial-translation invariance of the unforced system). 

As we are assuming that the reference state looses stability at fi = fiQ via a homoge- 
neous Hopf bifurcation in the absence of forcing, matrix Ai {^,k 2 ) must have a pair of 
complex-conjugate eigenvalues {Ai,A 2 } = {A (/x, k 2 ) , A (//, A; 2 )} (the overbar denotes com- 
plex conjugation) governing the instability, i.e.: 

(i) Close to the bifurcation ReAj> 3 < 0, while Re A can become positive for some k's, 

(ii) At the bifurcation Re A is maximum and null at k — (the perturbation with largest 
growth rate is spatially homogeneous): 

ReA = 0, (<9 fc ReA) = 0, (d 2 k Re\) Q < 0. (10) 

where, here and in the following, 

subscript affecting functions denotes {/i = /i , a = 0, k = 0} . 



(iii) The instability is oscillatory, i.e., 

Im Ao = ujo 7^ 0. 
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(11) 



Finally, the preceding properties imply that: 

(iv) All eigenvalues of .Mo = Jo, see Eq. 0, have negative real part but {Ao,Ao} = 
{iu , -iu }. 

For the sake of later use we introduce the right and left eigenvectors of Jo associated with 
eigenvalues {ioJo, —iuio}, 

Jo-h = iu h, J -h= -iu h, 
h f • Jo = iuoh^, ht • J = -iu W, 

where the short-hand notation h = Wi (/i = /io, k 2 = 0) , h = W2 (// = /io, k 2 = 0) has been 
introduced. These vectors verify the following orthonormality relations: 

h f - h = 0, h 1 " • h = 1, (13) 

as is triviall to be checked. 



IV. SCALES 

We are interested in determining the small amplitude solutions that emerge in the system 
close to the Hopf bifurcation, a parametric region that we define by 

fi = no + £ 2 /i2, (14) 

where e is a smallness parameter (0 < £ < 1). The study is based on the widely used 
technique of multiple scales J4J. These spatial and temporal scales appear naturally close 
to the bifurcation and are those on which the asymptotic dynamics of the unforced system 
naturally evolves. As is well known, in a homogeneous Hopf bifurcation these slow scales 
are given by 

T = e\X = ex, (15) 

which follow from the behaviour of A close to the bifurcation, Eq. (|14p . for values of k close 
to the most unstable mode k = 0: 

A (no + e 2 fi2, k 2 ) = A + (^A) e 2 ^ + \ {d 2 k \) Q k 2 + max {O (e 4 ) , O (k 2 ) } , (16) 
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where a term (<9fcA) k has not been included since (<9fcA) = as A is an even function of k. 
From this equation we obtain, making use of Eq. ([TO]) . 

Re A (/i + e 2 fi 2 , k 2 ) = (^ReA) e 2 /i2 - ~ \d 2 k Re A| k 2 + max {O (e 4 ) , O {k 2 ) } , 

what indicates that the only modes which can experience linear growth verify k = 0(e), 
hence the asymptotic dynamics of the system exhibits spatial variations on a scale x ~ 
k^ 1 ~ e^ 1 and the slow spatial scale X = ex follows. Thus, setting k = eki, 

1 



Re A (/x + £ 2 ^2, £ 2 kl) = e" 



(d,Re\) ^--\d 2 k Re\\k 2 l 



+ O (e 4 ) 



which shows that the growth (or decay) of the perturbations occurs on a scale t ~ (Re A)~ ~ 
e~ 2 and the slow timescale T = e 2 t follows. On the other hand, making use of Eq. (1111) and 
setting k = ek\ again, Eq. (fTBj) yields 



Im A (fi + e 2 fi 2 , e 2 kf) = u + e 2 



(d fM ImX) ^ 2 + - (d 2 k lm\) Q k 2 1 



whose first term, uq = O (e°), indicates that the original timescale t must be retained, while 
the rest of terms do not introduce other relevant timescales. 



A. Scales for the forcing 

As for the external forcing a (x, t) we assume that it is weak, of order e 2 , periodic in time 
at a frequency Uf close to the Hopf frequency lo , and depends on space on a scale of order 
£ -1 / 2 , which is shorter that the typical spatial scale X = ex of the system. We note that 
this choice introduces another spatial scale 

t = e l/2 x. (17) 

For the sake of definiteness we assume that U{ = u>o+e 2 u>2 and adopt the following expression 
for a: 

a (x, t) = e 2 (Xi (x, t) = e 2 A (£) exp {iu^T) exp (iuot) + c.c. (18) 

We note that the final result of the derivation does not depend on the harmonic character 
of the forcing: Any periodic (in fact almost periodic) a at frequency Uf close to ujq can be 
expressed as a Fourier series (with slowly varying coefficients) with fundamental frequency 
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a. 



Wo, of which (fl8|) is its first term, and the final result depends only on it [5|. We note for 
later use that 

F (/i, a) = F (^o + e 2 /i2, e 2 a 2 ) = e 2 a 2 (d a F) + O (e 4 ) , (19) 
as follows from ([5]). 



B. Scales for the system's oscillations 



Under all the previous conditions a multiple scale analysis is possible |4j and we look for 
asymptotic solutions to Eq. (0J in the form 

u (x, t) = V°° e m ' 2 M m/2 (X, T, £, t) , (20) 

^— — <m=2 

so that the expansion starts at order e, as usual. 

We finally introduce Eqs. ( Fl4|) and (TT8H20]) into Eq. (j4]) making use of the following chain 
rules for differentiation 



d t u = V°° e m ' 2 (d t + s 2 d T ) u m/2 , (21) 

<m=2 

V 2 u=V°° ^/ 2 (^| + 25 3 / 2 ^+5 2 9|)u m/2 , (22) 

^— <m=2 ^ 

and solve at increasing orders in e. Note that the occurrence of a term 2e 3 ^ 2 d^dx in the 
Laplacian operator suggests performing the expansion in terms of e 1 ^ 2 , and not of e, see 

(EDI. 



V. THE ASYMPTOTIC ANALYSIS 

The general form of Eq. (jlj) at any order e m is found to be 

3(u m ) = g m (X,T,£,t), (23) 

where 

3 (u m ) = <9 t u m - Jo ■ u m , (24) 

and g m does not depend on u m but on u n<m . Clearly, as 3 (exp (icjot) h) = 
Z (exp (— zwoO h) = 0, see Eq. (fT2"j) . the solvability of Eq. (j2"3"j) requires 

-t+27r/iJo 



2tt 



dt'ht • g m (X, T, £, t') exp (-zw f) = 0, (25) 
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(or its equivalent complex-conjugate) which ensures that g m does not contain secular terms 
(proportional to exp (iuot) h or to exp (—ioj t) h). Once condition f l25|) is verified, the asymp- 
totic solution to Eq. ( )23l) reads 

u m (X, T, f , t) = w m (X, T, f ) exp (iu; i) h + u m (X, T, f ) exp (-iu t) h 

+ u^(X ) T,£,t), (26) 

where u m (X, T, £) is not fixed at this order and the last term is the particular solution. Note 
that the solution f l26|) should involve, in principle, terms proportional to all the eigenvectors 
of 3 (•) [which are those of j7~o, see Eq. (jS}]. However all of them (but the first two) are 
damped according to exp [— |Re Aj (fjb , 0)| t], since Re Aj> 3 (ji = /x , fc = 0) < (and of order 
e° by hypothesis), except those associated with (h, h). 

A. Order e 

This is the lowest nontrivial order and 

gi = 0, (27) 

so the solvability condition ( 1251) at this order is automatically fulfilled. Then, according to 
Eq. (J26D, 

U! = ui (X, T, £) exp (iu t) h + u\ (X, T, £) exp (— ztu t) h, (28) 

where the scaled, slowly varying complex amplitude of oscillations u\ (X,T,£) is undeter- 
mined at this stage. We need to continue the analysis in order to meet solvability conditions 
that fix an equation for ui, which is our goal. 

B. Order e 3 / 2 

At this order we get 

g 3 /2 = 0, (29) 
hence the solvability condition is automatically fulfilled again and u 3 / 2 reads, according to 

U3/2 = U3/2 (X, T, £) exp (iu t) h + u 3/2 (X, T, ^) exp (-iu t) h. (30) 
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C. Order e 2 



Trivially one has 

g 2 = a 2 (d a F) + X? ■ <9f Ul + K (/i , 0; u 1; m) . (31) 

Making use of Eq. (1251) and taking into account that K is symmetric and bilinear in its two 
last arguments, Eq. ( 13T1) can be written as 



g 2 = « 2 (<9 Q F) + (X> • h) dim (X, T, exp (iwot) + (P • h) df«i (X, T, exp (-iu t) 

+ 2K (/x , 0; h, h) \ui\ 2 + K (/x , 0; h, h) w 2 exp (i2u;ot) + K (/j, , 0; h, h) w 2 exp (— i2cu t) . 

(32) 



The solvability condition (1251) reads now 

t-t+2ir/u)(i 

(ht • P • h) (X, T, + [ht • (0 a F) o ] / rft'« 2 exp (-i Wo O = 0, (33) 

which making use of ( 1T51) becomes 

c 2 9| Ul (X, T, + c 4 A (0 e^ T = 0, (34) 

where we defined 

c 2 = h f ■ X? ■ h, (35) 
c 4 = ht ■ (<9 a F) , (36) 

Condition (l34"j) implies that «i can be expressed as 

Ul (X, T, = (c 4 /c 2 ) MA (0 e^ T + ^ (X, T) , (37) 



where (£) is the particular solution to 

<*V (0 



^ (0 , (38) 



de 

and C/j (X, T) is a yet undetermined function of the slow scales X and T. Note that the 
solvability of Eq. ( 1351) requires that A does not contain a constant term. Hence we assume 
in the following that 

(A(O) = 0, (39) 
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where the angular brackets denote averaging over the spatial scale £. 
Finally, according to Eq. f l26|) . 



u 2 = u 2 (X, T, £) exp (iuj t) h + u 2 (X, T, £) exp (— iu t) h 
+ v |mi| 2 + v 2 u\ exp (2iu>ot) + v 2 u\ exp (— 2icu Q t) 
+ [v A A (0 exp (zwaT) + (Z>„ • h) df Ml (X, T, 0] exp (iw„*) 
+ [v A Z (0 exp (-iu 2 T) + (V ■ h) (X, T, 0] exp (-iw t) (40) 



where, 



v = -2J - 1 -K( / i ,0;h,h), (41) 
v 2 = - (Jo - «2o;oX) _1 • K (/i , 0; h, h) , (42) 
v A = - ( Jo - iujqI - iu h ® h f + iu h <g> ht j 1 ■ (<9 a F) (43) 

are constant vectors, and X is the NxN identity matrix. Note that both J§ and (j7o — i2ujqX) 
are invertible since neither nor 2iu$ are eigenvalues of Jq by hypothesis: otherwise other 
eigenvalues different from {iuo, — iuo} would have null real part at the bifurcation. For 
the same reason (^J — iloqI — iwoh ® ht + icooh eg) ht j is invertible too since iu is not an 
eigenvalue of (j — ico h <g> ht + iu h <g) ht j because we removed the subspaces spanned by 
{h,h}. 

D. Order e 5 / 2 

At this order we get 

g5/2 = 2<%<9 x ui + <9|u 3/2 + K (/i , 0; ui, u 3/2 ) = <9|u 3/2 + K (/i , 0; u 1; u 3/2 ) , (44) 

where the last equality comes from ( f28l) . The solvability condition ( |25l) reduces in this case 
to 

/t+2ir/u) 
dt'h^ ■ <9|u 3 / 2 exp (—iu t') = 0, 

as K (/x , 0; u 1; u 3 / 2 ) does not contain terms oscillating as exp (±ioj t). Making use of (!30|) 

we get 

<9|u 3/2 (X,T,O=0. (45) 
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Then, either u 3 / 2 = or w 3 / 2 does not depend on £. In the following we will consider the 
more general case, 

U3/2 = M3/2 (X, T) exp (iu t) h + u 3/2 (X, T) exp (-iu t) h. (46) 
Once the solvability condition has been verified, u 5 / 2 can be written, according to fl26|) . 

as 

U5/2 = M5/2 (X, T, f ) exp (zo; t) h + U5/2 (X, T, £) exp (-iu t) h 

+ u 5 J - /2 (X,T,£,t), (47) 

where the last term can be computed easily but we will not do as we need not knowing its 
expression. 

E. Order e 3 

This is the final order to be considered as it provides the shought equation for the small 
amplitude of oscillations. At this order we find 

g3 = -9 T ux + /i 2 (d^J) Q ■ ux 

+ V ■ d x ut + 2£> • dxd^u 3 / 2 + T> • <9f u 2 

+ 2K (/i , 0; u l5 u 2 ) + K (// 0) 0; u 3/2 , u 3/2 ) + L (/i , 0; uj, m, Ui) . (48) 

Note in the previous expression that dxd^u 3 / 2 = according to ( 14"6"|) . Application of the 
solvability condition ([25]) yields, after substituting Eqs. (125)1 . (jlB) and (HOD into Eq. (gSJ , 
making use of the symmetry and linearity properties of vectors K and L, and after simple 
but tedious algebra, 

d T ui = cifi 2 ui + c 2 d x ui + c 3 |mi| 2 u\ + c 2 d\u 2 + c A d\A (£) e iW2T + c 5 <9|wi, (49) 
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where 



Cl 




(d u J) n ■ h, 


(50a) 




— Vit 

— IT • 




^OUD ) 


c 3 


= 2h f 


• K (jiq, 0; h, v ) + 2h f • K (// , 0; h, v 2 ) + 3^ • L (/i , 0; h, h, h) , 


(50c) 


c 4 


= ht- 


(0«F) O , 


(50d) 


ca 


= ht- 


Do • va, 


(50e) 


c 5 


= ht- 


D • D ■ h, 


(50f) 



are constant coefficients. 

VI. THE COMPLEX GINZBURG-LANDAU EQUATION 

There remains just substituting ( 157)) into (|49|) . Once this is done, one can see that in Eq. 
( 149)) there are terms depending on the slow spatial scale X as well as terms depending on 
the "fast" spatial scale £ so that Eq. ( H9l) can be split into two equations: one containing 
just functions of X alone, and one containing functions of £ {and possibly X). The terms 
depending on £ determine partially the value of u 2 . The terms depending only on X deter- 
mine the equation of motion for Ui, see Eq. ( 157)) . which is the leading order amplitude of 
oscillations. 

The result is 

d T U x = d^Ur + c 2 d\U x + c 3 \U x f U x + c 3 (c 4 /c 2 ) 2 ^^Ux + 2c 3 \c 4 /c 2 \ 2 JU X , (51) 

7 = <^(o>, y=<i«A(or>, (52) 

where () denotes a spatial average (over the scale £) as already introduced. In order to arrive 
to this equation we assumed that {\ua {Of u a (£)) = 0. 

It is convenient to remove the explicit time dependence in (j5ip by performing the following 
change, 

U = U ie ~ iuJ2T } (53) 

and the equation becomes 

d T U = (ci M2 + zw 2 ) [/ + c 2 d 2 x U + c 3 \U\ 2 U + c 3 (c 4 /c 2 ) 2 7 F + 2c 3 |c 4 /c 2 | 2 7 'f/, (54) 
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which is of the type ([T]) with n = 2 as anticipated. 

Summarizing, for a system like (j2J) close to a homogeneous Hopf bifurcation, defined as 

fx = Ho + e 2 /i 2 , (55) 

and under the type of forcing analyzed along this paper, namely 

a (x, t) = e 2 A (£) exp (ia^T) exp (iu>ot) + c.c, (56) 

small oscillations emerge in the form 

(X, T, exp (iw t) h + e«i (X, T, exp (-iw t) h + O (e 3/2 ) , (57) 

where 

U! = [(c 4 /c 2 ) u A (£) + 17 (X, T)] e^ 2T , (58) 
ma (^) just follows passively the forcing via 

^ = -MQ, (59) 

where the forcing amplitude must verify (A (£)) = 0, and C/ is governed by (|54l . in which the 
different coefficients are defined in ( 1501) and in ( 1521) . Finally the validity of ( )54"|) requires that 
(I m a (0| 2m a (0) = 0- Note that all the conditions imposed on A (or on ua) imply that the 
sign of A should alternate in space. As an example, all previous conditions on the forcing 
are met for the simple case A(£) = A cos (<?£), in which case Ua{£) = (A Q /q 2 ) cos and 
7 = l' = \ (A /q 2 f. 

Finally, Eq. ( 154|) is valid whenever Rec3 < (supercritical bifurcation) since otherwise 
it would lead to unbounded solutions. If Re c 3 > the bifurcation is subcritical and the 
analysis must incorporate higher orders in the e-expansion. 



VII. A REMARK ON THE VALIDITY OF THE CLASSICAL COMPLEX 
GINZBURG-LANDAU EQUATION WITH RESONANT FORCING TO THE 
PRESENT CASE 

Equation (1541) is a CGLE with broken phase invariance, because of the occurrence of 
the term proportional to U. This term is typical of self-oscillatory systems forced at twice 
the natural frequency (2:1 resonance), as stated in the Introduction, but here it has been 
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obtained for a forcing resonant (1:1 resonant in fact) with the oscillations, whose amplitude 
is spatially modulated at a " short" spatial scale (shorter than the typical one in the unforced 
case). Why has this happened? 

Coming back to the 1:1 resonance its universal description is given by the CGLE ([I]) with 
n = 1, i.e. 

d t u = a\u + ci2V 2 u + a 3 \u\ 2 u + a 4 . (60) 
The derivation of this equation requires formally that 

a = e 3 a 3 = e 3 A 3 (X, T) e iuJot + c.c. (61) 

of order e 3 , not e 2 as we assumed up to here, and then 04 is proportional to A 3 (X, T), see 
j^]. The question is then: Is Eq. (1601) valid even when forcing is a bit stronger (of order e 2 ) 
and acts on shorter spatial scales, as we are considering in this paper? In order to make a 
closer analysis we consider the version of Eq. ( 1601) adapted to our notation, as derived in 

B, 

drUi = ci/i 2 «i + c 2 d 2 x ui + c 3 \ui\ 2 u x + c 4 A 3 (X, T) , (62) 

where all coefficients have the same meaning as in Eq. (1541) . 

Assume now that A 3 is "large" (this does not mean that the actual forcing a is) and that 
it varies on a "short" spatial scale £ = e~ 1 / 2 X, see f)17p and f)15p . In particular we consider 
formally that 

A 3 = e- x A (e-^X) e iuJ2T , (63) 

where the exponential e tuJ2T has been included in order to consider the same case we have 
been dealing with, namely that the frequency of forcing is tUf = uo + e 2 U2, see (I6T!) and (|T8j) . 
Note that with scaling (163]) we are formally in the same situation as in ffl8|) . What we will 
show next is that, if in Eq. (1621) we consider formally the scaling (1631) for the forcing, even 
if, apparently, this scaling is at odds with the valifity conditions applicable to that equation, 
one obtains the very same Eq. (1541) we have obtained before. This means that Eq. (I62p can 
be regarded as valid even when forcing is "strong" and varies on a "short" spatial scale. 

The derivation is reasonably simple: As a new spatial scale has been introduced we 
assume that Ui in fl62|) can be written as 

ui (X, T) = u? (X, T, + e 1/2 «S 1/2) (X, T, + e'u? (X, T, + O {e 3 ' 2 ) , (64) 
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what implies that the Laplacian will act as d\ — > d\ + 2e^ 1 / 2 dxd^ + e^d^. Then the 
asymptotic analysis starts. At the leading, e^ 1 order we obtain dju^ = — (C4/C2) A (£) e tul2T , 
what implies that 



it 



{0) (X, T, = (04/02) u A (0 e^ T + I7i (X, T) , (65) 



where 

d 2 u A (0 



-A (0 , (66) 



which are nothing but Eqs. (1371) and (1551) . respectively. 
The next, £ _1//2 , order reads 

= 2d x d 6 uf\ (67) 
which is identically fulfilled because of (1B"5]) . Then, at order e° we get 



arfJ = c lf i 2 u? + c 2 (d\uf + 2d x d^ /2) + dfuU) + c 3 J ^ 0) . (68) 
Substitution of (1631) leads to 

<9 T tA - CxMsC/i - c 2 9|^i - c 3 |f/i| 2 f/i - 2c 3 |c 4 /c 2 | 2 iU x - c 3 (c 4 /c 2 ) 2 7 e 2 ^ T Z7 1 
= c 3 [2 |c 4 /c 2 | 2 (\u A \ 2 - 7') ^ + (c 4 /c 2 ) 2 (t£ - 7) e 2 ^ r C7i] 
+ (C1M2 - iw 2 ) (c 4 /c 2 ) n A e^ 2T + c 2 (29 x 9^S 1/2) + dfu^) 

+ c 3 [|c 4 /c 2 | 2 (04/02) \u A \ 2 u A e^ T + (c 4 /c 2 ) ^e^t/ 2 + 2 (c 4 /c 2 ) M A e^ T l^ 2 ] , (69) 

where 7 and 7' have been defined as in ( l52l) . We note that in this equation the left hand 
side is independent of the short spatial scale £. On the other hand no term on the right 
hand side depends solely on X (we assume that (\u A \ 2 u A ) = as we have done in the rest 
of this paper). Then the solution to (169]) is simple: equate to zero both sides. Doing it with 
the right hand side gives a condition on 2d x d^Ui +d^u[ , which we are not interested in. 
Equating to zero the left hand side yields the sought equation, 

d T Ui = Cl fi 2 Ut + c 2 d 2 x V x + c 3 U x + c 3 (c 4 /c 2 ) 2 7 e 2 ^ T !7 1 + 2c 3 |c 4 /c 2 | 2 7^1, (70) 

which coincides exactly with (l5ip. 

This demonstrates that the usual CGLE describing the 1:1 resonant forcing of self- 
oscillatory systems, Eq. ( )62|) . is valid even when the forcing term is "large" and varies 
on a short spatial scale. The reason for this is quite easy to understand: Equation (1621) is 
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valid for positive \ii (above the bifurcation), for negative \l 2 (below the bifurcation), but 
even for fi 2 = 0, as is trivial to be checked. Then one should consider that Eq. f|62l) is 
valid whichever the value of \x-i be (it must be however, at most, of order e°). Then a trivial 
rescaling of §2$ with A 3 (X, T) = A (X) e^ T , 

t = V T, Z = V 1/2 X, ^^,r)=r ] - 1 / 2 u 1 (X,T), r / = /i 2 Rec 1 , (71) 

leads to 

d T ^j = (l + i6)^ + c 2 d 2 ^ + c 3 |^| 2 V + c 4 ?T 3/2 ^3 {v~ l/2 Z) e iri ~ lul2 \ (72) 

where 9 = Then, we see that there exists a normalization telling us that, if jj 2 is 

small (then rj is too) the forcing term in fl62|) can be effectively large and vary on a short 
spatial scale. In this case the detuning u 2 must be accordingly small, of order rj, so that 
f]~ 1 u>2 = O (f]°), as otherwise the inhomogeneous term is highly nonresonant and produces 
no effect. 

VIII. CONCLUSIONS 

Starting from a general, unspecified model for a spatially extended system affected by a 
homogeneous Hopf bifurcation and forced externally by a perturbation that is resonant in 
time with the system's oscillations (1:1 resonance) and is spatially modulated (involving sign 
alternations) at a scale shorter that the typical one of the unforced system, I have shown that 
the close to threshold dynamics of the system is governed by a complex Ginzburg-Landau 
equation with a phase symmetry breaking term, proportional to the complex conjugate of the 
oscillations amplitude. This term appears in the universal description of the 2:1 resonance 
of self-oscillatory systems (forcing at twice the system's naural frequency) and is responsible 
for the emergence of phase bistability and phase bistable patterns, including phase domain 
walls, bright solitonic structures and periodic patterns. Thus the kind of forcing put forward 
in this paper represents an alternative to the classical 2:1 periodic forcing. 

Finally I have shown that even the classical complex Ginzburg-Landau equation valid in 
a 1:1 resonance, which contains an inhomogeneous term (that can vary, in principle, only 
on long spatial and time scales), is valid to describe the phenomenon analyzed here, i.e., 
in that equation one can consider that the forcing is "large" and that varies on a "short" 
spatial scale. Even if not shown here, the same applies if the forcing term varies on a "short" 
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temporal scale, a case that has been analyzed in the context of "rocking" jfj], a 1:1 resonant 
forcing technique in which the amplitude of forcing is spatially uniform but varies in time, 
which has been studied theoretically and experimentally in different contexts {7 10] 
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